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In the layered transition metal oxide LiNb02 the Nb''"'' (4d^) ion is trigonal- prismatically coor- 
dinated with O ions, with the resulting crystal field leading to a single band system for low energy 
properties. A tight-binding representation shows that intraplanar second neighbor hopping t2 = 100 
meV dominates the first neighbor interaction {t\ = 64 meV). The first and third neighbor couplings 
are strongly modified by oxygen displacements of the symmetric Raman-active vibrational mode, 
and electron-phonon coupling to this motion may provide the coupling mechanism for superconduc- 
tivity in Li-deficient samples (where Tc « 5 K). We calculate all zone-center phonon modes, identify 
infrared (IR) and Raman active modes, and report LO-TO splitting of the IR modes. The Born 
effective charges for the metal ions are found to have considerable anisotropy reflecting the degree to 
which the ions participate in interlayer coupling and covalent bonding. Insight into the microscopic 
origin of the valence band density, composed of Nb d^2 states with some mixing of O 2p states, is 
obtained from examining Wannier functions for these bands. 

PACS numbers: 71.20.-b, 71.20.Nr, 71. 23. An, 74.70.Dd, 63.20.Dj 



I. INTRODUCTION 

The quasi-two-dimensional (2D) compound 
LiNb02 has become of interest in condensed matter 
studies in no small part due to the existence 
of superconductivity-'^ in the de-lithiated phase 
Li2:Nb02. This system consists of a triangular 
lattice of transition metal (niobium) ions separated 
by layers of O ions from 'blocking' (Li) layers, 
thus possessing structural similarities to high 
temperature superconducting cuprates. Important 
differences from the cuprates include the presence of 
a 4o? rather than 3d ion, and of course the triangular 
rather than square lattice. The trigonal prismatic 
coordination of Nb with six O atoms causes the 
valence band to be composed solely of Nb A^-^ states. 
Other d states are well separated by a 2 eV gap, 
leaving only a single band per formula unit required 
for analysis. The comparative structural simplicity 
of LiNb02 makes it attractive for theoretical stud- 
ies, especially as it requires modest computational 
time to characterize its properties. 

This alkaliniobate is also of interest due to the 
increasing attention garnered by niobates in recent 
yearsii^ The three-dimensional LiNbOa is an impor- 
tant and well studied ferroelectric material.'^ Super- 
conductivity at Tc ~ 20 K has been reported by a 
few groups^i^ in the barium niobate system. The 
phase has not been positively identified but it has 
been suggested to be a nonstoichiomctric BaNbO^; 
cubic perovskite. Superconductivity in the 5-6 K 
range has been reported in the Li-intercalated lay- 
ered perovskite system Li2;AB2Na„_3Nb„03„+i [A 
= K, Kb, Cs; B = Ca, Sr, Ba; n = 3 and 4]^^ 
This system has clear structural motifs in common 
with the peculiar Nbi2 029 system, where the 'single 
deficiency' (with respect to insulating Nb205 



= Nbi203o) releases two carriers into the cell; one 
localizes and is magnetic while the other is itinerant 
and results in conducting behavior^iiS 

While there is good agreement as to structure type 
and lattice parameters of the stoichiometric phase 
LiNb02, there is some disagreement in the literature 
as to its electronic and magnetic properties. Orig- 
inally it was reported as being paramagnetic with 
strong field dependence, however calculations pre- 
dict it to be a diamagnctic semiconductor.^^ Other 
experimental work suggests that it is diamagnetic, 
and strongly field-dependent paramagnetic proper- 
ties are the result of impuritiesJi^ The values of the 
band gap are also in mild dispute. Extended Hiickel 
calculations with a single isolated Nb02 layer give a 
bandgap value of 1.4 eVji^ and full potential linear 
muffin tin orbital calculations with the full crystal 
structure of LiNb02 give 1.5 eVJ^ Measured opti- 
cal reflectance shows an onset at ~ 2 eV, which is 
suggested to be a direct bandgap responsible for the 
burgundy-red color of LiNb02.^'^ Band structure cal- 
culations given here and by Novikov et ali^ are in 
agreement that LiNb02 is a direct gap semiconduc- 
tor. 

Geselbracht, Richardson and Stacji first reported 
that the de-lithiated phases Lia;Nb02 with x = 
0.45 and x — 0.50 support superconductivity at 
temperatures below 5.5 K and 5 K, respectively. 
Moshopoulou, Bordet and Capponii^ later reported 
samples with x ranging from 0.69 to 0.79 showing 
the onset of the Meissner effect at 5.5 K, but samples 
with x = 0.84 and above do not exhibit any super- 
conducting transition down to 2 K. Superconduc- 
tivity has also been reported in hydrogen-inserted 
H2:LiNb02, in which both samples reported had the 
same Tc — 5.5 K for x = 0.3 and x = 0.5i^ 

Removal of Li has the effect of adding holes to 



the conduction band made up of Nb d^^ states. As 
mentioned above, there does not appear to be any 
notable dependenc o^'^^ of Tc on the concentration 
of Li in the range 0.45 < x < 0.8. Such concentra- 
tion independence has also been observed in the 2D 
electron-doped system.^^-^* Lia;ZrNCl. For a phonon- 
paired superconductor with cylindrical Fermi sur- 
faces this behavior is understood in terms of phase 
space restrictions on phonon scattering processesiS, 
in 2D. The independence of Tc on carrier concen- 
tration in this system is quite surprising however, 
since the density of states (DOS) plots presented by 
Novikov et al}'^ and here varies strongly with energy, 
whereas parabolic bands lead to a constant DOS in 
2D. There has also been some work done investigat- 
ing the ion mobility of LiNb02 to assess its potential 
usefulness as a battery material, 

In this work we study the electronic structure 
and bonding of stoichiometric LiNb02 using den- 
sity functional theory in the local density approxi- 
mation (LDA) . Density functional perturbation the- 
ory is used to obtain zone center phonons and Born 
effective charges. The structure contains a single pa- 
rameter z that specifies the position of the O layers 
relative to the Li and Nb layers, and which seems 
to be of particular importance in understanding the 
electronic bonding and electron-phonon coupling in 
this compound. We investigate the corresponding 
Raman active vibrational mode (beating oscillation 
of the O planes) and the effect this has on the band 
structure. We obtain a tight-binding (TB) model 
that involves many neighbors and helps one to un- 
derstand how this variation affects interactions be- 
tween the Nb atoms. Many aspects of the electronic 
structure, such as chemical bonding and effective 
charges, can be seen to be interrelated by examining 
the Wannier functions of the valence bands. 



II. STRUCTURE 

Structural information of LiNb02 has been re- 
ported by several experimental groups, leading to 
the assignment of space group PQ^/mmc (No. 194), 
with Li occupying sites 2a (0,0,0) with symmetry 
3m, Nb occupying sites 2d (|, -j, j) with symmetry 
6m2, and O occupying sites 4f (i, |, z) with symme- 
try 3m. The O sites include an internal parameter z 
specifying their height relative to the Li layer along 
the c-axis. The Li sites are centers of inversion; the 
Nb sites have z-reflection symmetry in addition to 
the 3to triangular symmetry within the layer but do 
not lie at centers of inversion. 

There is good experimental agreement on lattice 
parameters, a = 2.90A, c = 10.46Aiii-iii^ The dis- 
tance a between Nb atoms is quite close to the near- 
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est neighbor distance of 2.86 A in elemental bcc 
Nb, suggesting that direct Nb-Nb coupling should 
be kept in mind. However, Nb"^+ Ad orbitals will 
be smaller than in neutral Nb. In later sections, 
we will discuss details of the Nb-0 interactions, and 
show that the electronic properties of LiNb02 are 
very sensitive to changes in the Nb-0 distance. The 
vertical distance between Nb and O planes can be 
calculated by (i — z)c, and the distance between Nb 

and O atoms is given by ^ + (i ~ z)^c'^. 

Some disagreement exists about the value of z. A 
larger value of z indicates the O layers are closer to 
the Nb layers. Meyer and Hoppe^"'^ found z = 0.1263, 
Moshopoulou, Bordet and Capponii^ report z — 
0.128, Geselbracht, Stacy, Garcia, Slibernagel and 
Kweili report z = 0.12478, Tyutyunnik, Zubkov, 
Pereliaev and Kar'kinSi report z = 0.1293. These 
variations amount to a difference in O layer position 
of about 0.05A. 

Optimization of unit cell geometry within LDA 
leads to a = 2.847A, c = 10.193 A, z = 0.1211. This 
differs from the measured value by Aa « —1.8%, 
Ac w —2.5%, for a total volume discrepancy of 
AV/V « —6%. The calculated overbinding of the 
unit cell is somewhat larger than is typical in LDA 
calculations, perhaps due to the low (practically 
zero) valence electron density in the Li layer. In 
this regard, we note that our Nb pseudopotential 
does not include the Ap states in the valence bands, 
which may have some effect on structural proper- 
ties. Optimization with respect to z, holding the 
lattice parameters a and c fixed at experimental val- 
ues gives z = 0.125, which is very close to the value 
of Geselbracht et ali^ 



III. ELECTRONIC STRUCTURE 
CALCULATIONS 

The present results have been obtained through 
the use of the ABINIT codeM The lattice con- 
stants and the parameter z given by Meyer and 
Hoppeii were used in the calculations in this section. 
Pseudopotentials used were LDA, Troullier-Martins 
type^ generated using the Perdew-Wang LDA 
exchange-correlation functionali24 A self-consistent 
density was generated with 36 irreducible k-points, 
which was then used to calculate energies at 116 k- 
points along high symmetry directions to plot the 
band structure. The kinetic energy cutoff for the 
planewave basis was set to 40 Hartrees. 

Figure^of the important region of the band struc- 
ture shows the top of the O 2p bands, a 2 eV gap 
to two Nb Ad bands which are the highest occu- 
pied bands, and another 2 eV gap to the conduc- 
tion bands. The band gap, direct at the F point, is 




FIG. 1: (Color online) Band structure for LiNb02 with 
O height z — 0.1263 typical of the slightly scattered ex- 
perimental results. Two Nb bands (one per Nb layer) lie 
within a 5.5 eV gap between the O 2p bands below and 
the remaining Nb 4d bands above. The tight binding fit 
of Tabled is shown with a dashed line. 



1.9 eV, somewhat higher than the 1.5 eV reported 
by Novikov et al. (also at the F point) using the full 
potential linear mufhn-tin orbital method (FLMTO) 
with the same structural parameters!^ LDA calcu- 
lations often underestimate band gaps, so we con- 
clude that it is most likely the experimental band 
gap is somewhat larger than 2 eV. Measured opti- 
cal absorption suggests that the band gap should be 
around 2 eV,-''« which is consistent with the reported 
dark reddish color. ^^'^'^ 

Projected density of states (DOS) calculations, 
presented in Fig. |21 were done using the tetra- 
hedron method projecting to atomic spheres with 
radii of R(Li) = 1.7 a.u., R(Nb) R(0) = 2.0 
a.u. These calculations confirm that the upper- 
most valence bands are largely Nb 4d character, with 
roughly 15% contribution from O 2p states, and vir- 
tually none from Li states. These projected DOSs 
are consistent with those calculated by Novikov et 
ali^ The trigonal crystal field splitting of the Nb 
3d states results in the orbital energy centered 4 
eV below the other Nb 4d states (two Cg doublets, 
based on {xz, yz} and {x^ — y^, xy}). The unit cell 
contains two symmetry-related Nb layers which re- 
sult in the presence of two Nb dz^ bands in the band 
plot. The splitting of the bands along F-M-K-F re- 
flects the strength of interaction between layers, as 
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FIG. 2: (Color online) Density of states for Nb and O 
atoms. For LiNb02, states are occupied up to energy 
E = 0. The dominant states in the valence band for Nb 
are 4d-orbitals, and for O they are 2p-orbitals. Dashed 
vertical lines are included to indicate the Fermi energy 
in the rigid band picture for Lia;Nb02 for concentrations 
X = (half-filled band), x = 0.4, x = 0.8 from left 
to right. The last two values of x roughly bracket the 
concentrations where superconductivity is observed. 



discussed below. 

As a result of the large trigonal crystal field split- 
ting, the Nb dz2 state forms a single band, trian- 
gular lattice system of the type that is attracting 
renewed interest. The calculated DOS in Fig. [3 has 
only a vague resemblance to the nearest neighbor 
TB model (see for example Honercamp™) making it 
evident that longer range hopping must be impor- 
tant. We address this question in detail in section 

EH 

This general electronic structure, when hole- 
doped, is rather remarkable. It is a bona fide rep- 
resentation of a single band system on a triangular 
lattice, moreover the band is well isolated from other 
bands. Although single band models are widely used 
to investigate concepts and to model many-body ef- 
fects, there are very few actual realizations in real 
solids. The value of the effective Coulomb on-site 
repulsion U for Nb'^^ is not established, but with 
W < 1.75 eV it will not require a very large value to 
lead to correlated behavior. The hole-doped system 
Lia;Nb02 provides a unique platform for the detailed 
comparison of single band models with experimental 
data, which for thermodynamic properties and for 
spectroscopic data below 2 eV should involve only 
the single active band. 




FIG. 3: (Color online) Isosurface plot of the F-pt wavefunctions for LiNb02 for (a) the lower of the pair of valence 
bands, and (b) the upper valence band, made with XCrysDeni— Both wavefunctions are generated with the same 
isosurface value. Nb atoms are turquoise (grey), O atoms are dark red (dark gray), Li are light gray. Yellow and blue 
(light and dark) indicate opposite signs of the wavefunction. The symmetry around the Nb atom is apparent as 
well as the contribution to the bonding made by the O atoms. Trigonal prisms with O atoms at the corners and Nb 
atoms in the center have been outlined as a guide to the eye. 



IV. ZONE CENTER WAVE FUNCTIONS 

We can examine Kohn-Sham wavefunctions to 
gain insight into the electronic structure. Shown in 
Fig. Oare isosurface plots of the F wavefunctions for 
the lower and upper valence bands. Both plots reveal 
the (^22 -character of the bands on the Nb atoms, as 
well as 2pz bonding contributions from the O atoms. 
For the upper band, the p-lobes have more weight 
away from the sandwiched Nb layer, indicating the 
weight of bonding with O is somewhat decreased. 
The change in sign between the two Nb layers in 
Fig. I^Id indicates the presence of a nodal plane in the 
Li layer arising from the wavefunction's antibonding 
character (with respect to interlayer coupling). 

The presence of an interesting triangular shape of 
the wavefunction around the waist of the Nb atoms 
is more unexpected, being noticeably larger for the 
bonding function. The corners of these triangles 
do not point toward nearest neighbors, instead they 
point toward interstitial holes in the Nb layer which 
are not directly sandwiched between O atoms. This 
shape reflects either three-center Nb-Nb bonding, 
which would result in an increased density in this 
region, or antibonding character of direct Nb-Nb in- 
teractions that would decrease the density along the 
directions connecting Nb atoms. The character and 



origin of this triangular 'waist' around Nb atoms are 
clarified in the following section. 

V. WANNIER FUNCTION 

A single isolated band also provides an unusually 
clean system for performing the transformation from 
reciprocal space (bands) to real space (bonds) to 
look at character and strength of interatomic inter- 
actions from a local viewpoint. The straightforward 
way of doing this is to generate Wannier functions 
(WFs) defined by 

|Rm) = N^'^-" J2 C/,Le-*-^|kn) (1) 

where can be any unitary matrix, and A'k is the 
number of k-points used in the summation. The 
orthonormal Bloch states are denoted by |kn). If 
there were only one Nb02 layer per cell, hence only 
one band, U(k) would be simply a complex num- 
ber of unit modulus. Here it is a 2 x 2 matrix in 
the band space, which gives rise to very modest ex- 
tra complexity. These WFs display orthonormality 
(Rm|R'm') = 6u.,ii'Sm,m' ■ Once a choice for is 
made, this approach generates a basis set of identi- 
cal WFs for each point in the Bravais lattice of the 




FIG. 4: (Color online) Isosurface plot of the Wannier 
function for LiNb02, made with XCrysDenp— Niobium 
atoms are turquoise (gray), O atoms are dark red (dark 
gray), Li not pictured. Yellow and blue indicate opposite 
signs of the WF. The symmetry around the central 
atom is apparent, as well as the xy/x^ — character 
around nearest neighbors, which maintains the orthogo- 
nality of WFs centered on different atoms. 

system being studied. 

We have chosen the unitary matrix U such that 
U = M{M'^M)-^/^ with Mmn = (kn|5™) where 
\Sm) are a set of atom-centered functions with the 
desired c?22 symmetry. This technique has been used 
previously in the study of magnetic ordering in 
cupratesiSI Choosing \Sm) as hydrogenic or- 
bitals centered on the Nb atoms for the two entan- 
gled valence bands gives a pair of identical WFs for 
each Nb bilayer in the unit cell, based on the Nb 
cR, symmetry that was also observed in the F-point 
wavefunction. An 8x8x4 mesh totaling 256 k- 
points was used to generate the WF shown in Fig. 
21 Increasing the size of the mesh does not change 
the visual appearance of the WF. 

The WF naturally displays the prominent lobe 
along the ±z axis that is characteristic of the iz^ — r^ 
function. It also contains Pa contributions from the 
nearest O ions, as well as smaller but still clear 2p 
amplitude on the second neighbor O ions that is ap- 




FIG. 5: (Color online) Isosurface plot of the difference 
|W^z=o.i359(r)| - lWz=o.i263(r)|, looking along the [010] 
cartesian direction. Yellow (light gray) is the positive iso- 
surface and blue (dark gray) is negative, with Nb atoms 
in turquoise (gray) and O atoms in dark red (dark gray) . 
As the O atoms are moved closer to the Nb layer, there is 
an increase in density above and below the Nb atom; in 
addition the contributions from the O-p orbitals appear 
to tilt more towards the central Nb atom. 



proximately Ptt- Here Pa denotes the combination 
Px,Py,Pz functions oriented toward the reference 
Nb atom, while is perpendicular. The p^ lobe 
nearest the Nb ion has the same sign as the 3z^ — 
lobes along the z axis. 

The in-plane part of the WF reveals more about 
the bonding in LiNb02. The ring-structure of the 
— function in the x ~ y plane connects with 
combinations oixy and x^ — y^ functions on the three 
neighboring Nb ions to form a substantial in-plane 
windmill structure with three paddles, which are ori- 
ented between the ^-projections of the neighboring 
O ions. A "3^2 - r^" symmetry WF has the full 
symmetry of the Nb site, which allows a contribu- 
tion of an appropriate combination of x'^ — y^ and 
xy atomic orbitals, and the windmill structure re- 
flects the bonding of such functions on neighboring 
Nb ions. Thus the 4d character is not entirely d^i^ 
which clarifies the threefold symmetric shape around 
the Nb waist in the k = wavefunction discussed in 
the previous section. 

In addition, there is a 'hotdog' structure of op- 
posite sign connecting pairs of nearest-neighbor Nb 
ions. Of the two types of Nb-Nb pairs, these hot- 
dogs lie nearest to the first neighbor O ions. Small 
nearest-Nb "backbonds" in the x — y plane of both 
the windmill arms and the hotdogs are also evident. 



As we note below, the WF contains information 
about the Nb-Nb hopping parameters. The more 
distant parts of the WF function are a direct reflec- 
tion of the existence of such longer range interaction. 

This Nb Adz2 WF is similar to that for Ta 5dz2 in 
isostructural TaSe2i^^' A difference is that the near- 
est neighbor O 2p orbitals contribute more signif- 
icantly to the bonding in LiNb02, giving a larger 
nearest neighbor hopping. The lobes above and 
below the transition metal appear to be somewhat 
more pronounced in TaSe2 . We also observe changes 
toward such appearance in this Nb WF as the O 
layer height (2) is increased, forcing the O atoms 
closer to the Nb layers and increasing ti. Given 
in Fig. 13 is a plot of the difference in the magni- 
tudes of the WFs taken for the two different z val- 
ues of 0.1359 and 0.1263 (O displacements of roughly 
0.1 A.) There is an increase in the density of the d^^ 
lobes above and below the transition metal as well 
as a tilting of the O orbitals toward the central site, 
as the Nb-0 distance is reduced. Correlating this 
with the tight binding results from the next section 
suggests there is an increase in the role the O atoms 
play in the hopping to nearest neighbors as the Nb- 
O distance is decreased. The opposite effect occurs 
when z — 0.1167; the O orbitals tilt to become ori- 
ented more along the z-direction, indicating the O 
atoms are interacting less strongly with Nb. 



VI. TIGHT BINDING 

The valence bands seen in Fig. ^ are well iso- 
lated from bands above and below. These bands 
contain the superconducting electrons in de-lithiated 
Li2:Nb02, so a simplified tight binding fit using only 
these bands would be useful in doing studies on sim- 
ple models. The TB model to fit two bands corre- 
sponding to identical Nb d^^-O p hybridized Wan- 
nier functions on the two Nb atoms per cell gives 
dispersion relations from the roots of 



(2) 



where e^^ arises from the hopping processes within a 
plane of atoms and from hopping between planes. 
Within the plane of the triangular lattice, the near- 
est neighbor dispersion for the six nearest neighbors 
is 
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(3) 

We have found that a good fit to the Nb bands 
requires several neighbors. Third and fifth neigh- 




FIG. 6; (Color online) Band structure for O layer po- 
sition z = 0.1167, a larger Nb-O layer separation by 
0.1 A compared to the bands in Fig. The band width 
is 25% smaller, and band gap is significantly less, com- 
pared to the experimental structure (z = 0.1263). 



bors in the plane can be included with this form by 
noticing odd neighbors just require a larger lattice 
constant. Second neighbors require a larger lattice 
constant and a 30° rotation. There are twelve fourth 
neighbors, which can be treated as two lattices with 
separate rotations. 

First and second neighbors in adjacent layers can 
be included with terms of the form 



4 = 2cos^^2t^ 
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It is worth noting that on the zone surface (where 
kz = ±f ), all terms in e^, including further neigh- 
bors neglected from Eq. Q are identically zero, re- 
sulting in exact degeneracy of the two valence bands 
at these zone surfaces. 

Using the Wannier basis given in the previous sec- 
tion, the hopping parameters can be directly calcu- 
lated (rather than fit) as matrix elements of the LDA 
Hamiltonian: 
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TABLE I: Tight binding hopping parameters for three values of the O internal structural parameter z. Hopping 
parameters are reported in meV; bandwidths and band gaps are reported in eV. Distance between Nb and O atoms 
are reported in angstroms. The largest hopping coefficients are in boldface. Squeezing the O layers closer to the Nb 
layers {z = 0.1359) results in a dramatic increase in ti, making it larger than t2, and a large change in as well. For 
z = 0.1359, the band gap is indirect. The direct gap at F is listed followed by the indirect gap in parentheses. See 
figure |7| is slightly refined, to make the splitting at F agree with the DFT result better. 











V 

















r M K FA L H A 



FIG. 7: (Color online) Band structure for O layer posi- 
tion z = 0.1359, about 0.1 A closer to the Nb layer than 
in experiment (Fig. 0. The bandwidth is 15% larger 
than for the experimental structure, and the peaks at 
M and K in the valence band are similar to the band 
structure calculation of Lio.5Nb02 shown by Novikov et 
al.^ suggesting that the underlying electronic structure 
of Lia;Nb02 is determined by the position of O layers 
rather than the Li concentration. 

Using the definition of the WF along with orthog- 
onahty of the eigenfunctions, Eqs. Q reduce to 

^|fl|,™ = A^k'E=k„|f7^„Pe'*'-^ (6a) 

kn 

^\R\,m,m' — -^k E]£kn(C^m'n) t^mn^''' ^- (6b) 
kn 

where R is the difference in centers of the WFs. De- 
viation of a perfect representation of the bands arises 



only from approximation of the summation in Eq. 
(6) or neglect of further neighbors. 

The results of the tight-binding calculation are 
presented in TableHJ as calculated from a 12 x 12 x 6 
Monkhorst-Pack grid of k-points for various values 
of z, along with values for 2H-TaSe2 as given in Ref. 
l28i Remarkably, the second neighbor hopping inte- 
gral t2 is the dominant contribution for z < 0.1263. 
Also notable is the strong dependence of first and 
third neighbor hoppings, ti and t^, on the Nb-0 
distance, which causes ti to become largest for z 
somewhat greater than 0.1263. The other param- 
eters in the table are all relatively independent of 
the Nb-0 distance. The third and beyond inter- 
plane hoppings may appear small, but they must be 
compared keeping in mind that there are twice as 
many third and fourth interplane neighbors as first 
and second. Since ti ^ 2t^ ~ 2t^, even these long 
distance hoppings are relevant. Indeed, it is not pos- 
sible to get good agreement in the band splittings at 
the points F, M, K, and the minimum between K 
and F without five interplane hopping parameters 
(four for z = 0.1167). 

The TB model for z = 0.1167 is strongly domi- 
nated by the second neighbor interactions between 
the Nb atoms, with other hoppings roughly an order 
of magnitude smaller. If t2 were the only hopping 
parameter present, the system would consist of three 
disjoint sub-lattices and thus states would be three- 
fold degenerate (although it would not be obvious 
on the standard band plot). Including ti as a per- 
turbation would couple them weakly and remove de- 
generacies. In this case, the unit cell of the system 
could be taken as the second-neighbor lattice and 
the Brillouin zone would be folded back so that the 
K(H) point in the band plots shown is mapped onto 
the F(A) point. The result of this folding back (and 
the splitting that occurs because of non-zero ti), is 
evident in Fig. El in that the eigenenergy of the H 
point is nearly the same as for the A point. The 
M(L) point is mapped onto itself, and the minimum 
between F(A) and M(L) is mapped onto the K(H) 
point, which results in a band structure that looks 



very much like a \/3 larger triangular lattice with 
only nearest neighbor hopping. Similar dominance 
of 2nd neighbor hopping has been observed for 2H- 
TaSe2 Those authors relate their small value for ti 
to phase cancellation of Wannier functions centered 
at first nearest neighbors, whereas in-phase contri- 
butions at second-neighbors give a large t2 matrix 
element. 

For z = 0.1359, we find that the nearest neighbor 
interaction becomes dominant. From Table and 
Fig. |S1 we see that the increase in hybridization of 
the O anion p-state in LiNb02 as the Nb-0 distance 
is decreased plays a significant role in first neighbor 
hopping, while at the same time this hybridization 
is rather insignificant in second neighbor hopping. 
This behavior suggests that a significant contribu- 
tion to the hopping integrals comes from within the 
Nb planes, perhaps through the hotdog structure 
noted in the previous section. The band structure in 
this case looks similar to that of Lio.5Nb02, calcu- 
lated by Novikov et al.j^'^ especially in the magnitude 
of the variations around the M and K points, relative 
to the bandwidth. Their result shows some splitting 
of the bands occurring at the zone surface (along 
the lines A-L-H-A), which may arise from the two 
Nb sites becoming inequivalent due to the removal 
of one Li atom. If the Nb sites were kept equivalent 
as the Li is removed, such as in the virtual crystal ap- 
proximation, then the similarities between the band 
structures suggest that the main impact on the elec- 
tronic structure is captured by changing the Fermi 
energy and allowing the O layers to move closer to 
the Nb layers, decreasing the Nb-0 distance. 

Another point of interest we see in Figs. and 
[3 is a significant change in how the valence band is 
situated in relation to the O 2p bands. Nb 4c? states 
in the conduction band around the K point are sig- 
nificantly lowered in energy as the Nb-0 distance is 
decreased so that the bandgap becomes indirect and 
much smaller. This change will have contributions 
from both a Madelung shift and from changes arising 
from the altered Nb-0 interaction. 



VII. ZONE CENTER VIBRATIONAL 
MODES 
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of SMODES, part of the Isotropy package of Stokes 
and Boyer.^^ Masses of the ions used were 6.94 (Li), 
92.91 (Nb), and 16.00 (O), aU in amu. Also in Table 
njare given values of the isotope shift a, defined by 



where i runs over the different (types of) atoms in 
the crystal. The derivative was approximated with 
a centered difference, which, for a harmonic crystal 
with uj oc Af^^/^, introduces errors of order (^)^- 
Using SM = 1 amu, this amounts to an error of 
about 2% for Li, and less for heavier elements. One 
would expect Ui = ^ for a harmonic mode. 

The dielectric tensor for a hexagonal system such 
as this will be diagonal with two distinct elements. 
For frequencies much higher than phonon frequen- 
cies but well below the gap, the static dielectric con- 
stants, calculated without ionic contributions, are 
found to be — 10.3 and — 3.78. One can then 
use the generalized Lyddane-Sachs- Teller relation, 
^0 = ^'^ nm«o,m/'^TO,m)^ (^here a runs over the 
cartesian directions, and the product is taken over 
all IR active modes with polarization in that direc- 
tion) to include the ionic contribution for the low 
frequency dielectric constant. Doing so, we find for 
LiNb02 that eg = 15.9 and eg = 8.27. It is typi- 
cal for LDA calculations to overestimate somewhat 
the dielectric constants, presumably due to the band 
gap being underestimated. 

The eigenmodes fall into groups which are eas- 
ily distinguishable by frequency. In the following 
discussion we include degeneracy in the numbering 
of modes. The first group, modes 1-3, are low fre- 
quency modes involving displacements of massive 
Nb-0 layers against each other, with Li layers re- 
maining fixed. The second group, modes 4-7, involve 
the lighter Li layers, and lie around 29 meV. The 
third group, modes 8-17 in the 53-65 meV range, in- 
volve motions where layers of atoms slide in the x—y 
plane against adjacent layers of atoms, or groups of 
layers bounce against each other in the z direction. 
The last group of modes at 85-95 meV involves lay- 
ers nearby each other bouncing against each other 
in the z direction. A more detailed discussion of the 
individual modes is given in Appendix 1X1 



We have calculated the 3N — 3 = 21 optical 
phonon eigenmodes with q = 0. To our knowl- 
edge, experimental IR or Raman measurements have 
not been reported on LiNb02. However, we be- 
lieve these results will be useful for comparison 
with electron-phonon calculations done on the de- 
lithiated Li2;Nb02, which will be explored in fu- 
ture work. The frequencies are given in Table ^ 
Mode symmetries have been obtained with the use 



VIII. EFFECTIVE CHARGES 

Born effective charge calculations have been car- 
ried out as described by Gonze and LeeM^ The cal- 
culation was done using a planewave cutoff of 60 Ha 
and 200 k-points in the Brillouin zone to give effec- 
tive charge tensors Z* for each atom which are diag- 
onal, and obey Z*^ = Z* , as required by symmetry 
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°IR active mode 
''Raman active mode 

TABLE XL Calculated zone-center phonon frequencies, with LO modes listed where they occur. See Appendix 1X1 for 
detailed descriptions of the atomic motions. All phonons with x-y polarization are doubly-degenerate, as required by 
the hexagonal symmetry. The last six columns show the calculated value of a^, defined by equation [7| for both TO 
and LO modes. 
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Li Nb 


Na Co 


7"* 

^zz 
^Formal 


1.10 2.26 -1.68 
1.69 1.31 -1.50 

+1 +3 -2 


0.87 2.49 -1.68 
1.37 0.87 -1.12 
+1 +3 -2 



TABLE IIP Calculated Born effective charges for 
LiNb02, compared with those calculated for NaCo02 
by Li, Yang, Hou and Zhu^ Note the disparity in the 
effective charges for the alkali metals and the transition 
metal in the z-direction. For O however, Z* is not far 
from isotropic. 



for a hexagonal crystal. Results are given in Table 
mil for the relaxed atomic structure. The effective 
charges of NaCo02 have been included in this table 
to illustrate the strong similarity between these two 
systems, which are similarly layered transition metal 
oxides but otherwise their electronic structures are 
quite different. 

The effective charges of LiNb02 are rather differ- 
ent than those reported by Veiten and Ghosez'^ for 
LiNbOs. There are structural and chemical differ- 
ences responsible for this; the structure of LiNbOs 
has Nb coordinated in a distorted octahedral envi- 
ronment and thus has low symmetry. Perhaps more 
importantly, the formal charges of Nb are very differ- 
ent: Nb3+ in LiNb02, but in LiNbOa, leaving 
the formal electronic configuration of Nb in LiNbOs 
as d'^. The effective charges of such ions are known to 
sometimes be rather large, as is the case in LiNbOa. 

Z*^(Li) is close to the formal charge of Li indi- 
cating ionic type response for in-plane motion. In 



fact, all ions have effective charge closer to the for- 
mal value for in-plane displacement than for z dis- 
placement. The charge tensor for Li shows similar 
anisotropy to that of Li in LiBC)^ where it was in- 
ferred that Li is involved chemically in coupling be- 
tween B-C layers. Similarly here, it can be expected 
that Li has chemical involvement in interlayer cou- 
pling between the electron rich Nb-0 layers more so 
than Na does in Na2;Co02. This involvement is not 
strong enough to inhibit de-intercalation, however. 

The smaller than formal charges for Nb indi- 
cates substantial covalent character to the bond- 
ing, which is especially strong in the z direction. 
This charge 'renormalization' may be related to the 
strong change in band structure due to change of the 
distance between the Nb and O layers, discussed in 
earlier sections. 



IX. SUMMARY 

In this work we have examined the electronic 
structure, the Wannier functions of the valence 
bands and their tight-binding parametrization, and 
several vibrational properties of the quasi-two- 
dimensional material LiNb02. These basic prop- 
erties are important in providing the basis for un- 
derstanding the behavior of the system when it is 
hole-doped by de-intercalation of Li. 

Li2:Nb02 seems to be a candidate for an insulator- 
to-superconductor transition, as it undergoes an 
insulator-to-metal transition and becomes super- 
conducting, but doping studies are not systematic 



enough yet to shed hght on this possibihty. Phase 
space considerations in two-dimensional materials 
lead to superconductivity that is independent of 
doping level, at least when the Fermi level is not 
far from a band edgeji^ This property supports the 
possibility of an insulator-superconductor transition 
with doping, i.e. as soon as it becomes metallic it is 
superconducting. 

The Wannier function, centered on and based on 
the occupied Nb d^^ orbital, illustrates graphically 
not only how Nb bonds with neighboring O ions, 
but also that it couples with oxygens in adjoining 
layers. Inspection of the k-space wavefunction sug- 
gests there are direct Nb-Nb interactions within the 
band as well as indirect coupling through O ions. 

A tight-binding representation of the band, ob- 
tained by direct calculation (rather than fitting) us- 
ing the Wannier function, reveals that several neigh- 
bors both in-plane and inter-plane are necessary to 
reproduce the dispersion. To model the interlayer 
interaction precisely requires four to five interlayer 
hopping parameters. For the experimental structure 
intralayer hopping is dominated by 2nd neighbor 
hopping, but the strength of both 1st and 3rd neigh- 
bor hopping is strongly modulated by varying the 
Nb-0 distance. This electron-lattice coupling pro- 
vides one possibility for the superconducting pairing 
mechanism in Lia;Nb02. 

Comparison with previous calculations of the 
band structure of Lio.5Nb02 suggests that the elec- 
tronic structure of Lia;Nb02 can be modeled in vir- 
tual crystal fashion, or perhaps even in rigid band 
fashion if that is convenient. The Nb d^^ band 
is remarkably well isolated from neighboring bands 
by Madelung shifts and crystal field splitting. We 
suspect that the overriding importance of this sys- 
tem will result from its unique standing as a single 
band, triangular lattice system for which sophisti- 
cated studies on simple models may be directly com- 
parable to experimental data. 



X. ACKNOWLEDGMENTS 

We have benefited from communication with K. 
E. Andersen and A. Simon. Some of the calcu- 
lations of LiNb02 wavefunctions used to compute 
tight-binding parameters were done on machines 
at Lawrence Livermore National Laboratory. This 
work was supported by National Science Foundation 
Grant DMR-0421810. 



10 

APPENDIX A: PHONON MODE 
DESCRIPTIONS 

Following is a descriptive list of the atomic mo- 
tions for each mode given in table ^ grouped by 
their spectroscopic activity. Degeneracy is included 
in the numbering of modes. 

1. Silent Modes 

Modes 1-2: Nb02 layers slide against each other 
in the x-y plane. 

Mode 3: Nb02 layers bounce against each other 
in the z direction. 

Modes 4-5: Li02 layers slide against each other 
in the x-y direction. Li displacements are roughly 
an order of magnitude larger than O displacements, 
which results in the mass of the O atoms having very 
little effect on the frequency of this mode. 

Modes 11-12: Adjacent O layers slide against each 
other, with Li and Nb layers remaining stationary. 

Mode 15: Li02 layers bounce against each other 
in the z-direction. This is the one high frequency 
mode whose isotope shift a gets its major contribu- 
tion from the Li atoms. 

Mode 20: Out-of-phase Li layers beat against ad- 
jacent O layers in the z direction. 

Mode 21: Out-of-phase Nb layers beat against ad- 
jacent O layers in the z direction. 



2. IR Active Modes 

Mode 6-7: In-phase Li layers sliding against Nb02 
layers. The main contribution to a comes from the 
light Li atoms, with the heavy Nb02 layers making 
very little contribution. This mode is also IR active, 
with a small LO-TO splitting Aw = wlo~'^to = 3.6 
meV. 

Mode 8: The TO mode consists of Li layers bounc- 
ing against Nb02 layers in the z direction, and can 
be thought of as the "z-polarized version" of modes 
6 and 7. When the macroscopic electric field is in- 
cluded, the O displacements change sign and the 
mode becomes Li-0 layers bouncing against Nb lay- 
ers. Acj = 16.7 meV. 

Modes 13-14: Li and Nb are in synchronized slid- 
ing against O layers. The metal layers are all in 
phase with each other, and out of phase with the O 
layers. This mode is IR active, with Acj = 7.2 meV. 

Mode 19: The transverse mode consists of Li02 
layers beating against Nb layers in the z-direction. 
The Li-0 layers are all in phase. This mode is IR 
active, and once the macroscopic electric field is in- 
cluded, this mode can be characterized as having Nb 



layers bounce against adjacent O layers, where the 
Nb layers are out of phase. Auj = 10.8 meV. 

3. Raman Active Modes 

Modes 9-10: Out-of-phase Li layers slide in the x-y 
direction against adjacent O layers. The rest of the 
phonon modes from here on have the contribution 
to their isotope shift a dominated by O. 

Modes 16-17: Out-of-phase Nb layers slide against 
adjacent O layers in the x-y plane. 

Mode 18: The O layers beat against each other 
in the z-direction. This fully symmetricRaman ac- 
tive mode is of particular interest, as it corresponds 
to the variation of the internal structural param- 
eter z. This mode has a significant impact on the 
electronic structure as discussed in previous sections, 
so we suspect that it may have particularly strong 
electron-phonon coupling in hole-doped Lia;Nb02. 
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the macroscopic electric field. Without the electric 
field, the eigenmodes can be described as one where 
Nb layers bounce against Li-0 units and one where 
Li layers bounce against Nb-0 units. The applica- 
tion of the electric field causes these vectors to mix, 
and the descriptions change so that the latter one is 
replaced by a mode where Li and Nb bounce against 
O layers. The distinction between these two modes is 
the relative sign for Li and Nb displacements (neg- 
ative for the former, positive for the latter.) It is 
not unusual to have modes mix to give different LO 
modes, as the non-analytic part of the dynamical 
matrix which incorporates the electric field can have 
eigenvectors differing from the analytic part of the 
dynamical matrix, resulting in a mixing of modes 
of the same symmetry. This is discussed in some 
detail by Gonze and Lee.'^^ This complicates the as- 
sociation of LO modes with their corresponding TO 
mode, so we have chosen to pair the LO and TO 
modes such that wlo > wto- 



4. Clarifying Comments 

The z-polarized A2„ IR active modes 8 and 19 
have their eigenvectors affected by the inclusion of 
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